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Abstract. We extend earlier calculations of the attenuation suffered by 7-rays during 
their propagation from extragalactic sources, obtaining new extinction curves for 7-rays 
down to 10 GeV in energy, from sources up to a redshift of z = 3. 

The recognition that high energy 7-rays, propagating over cosmological distances, 
suffer electron-pair-producing interactions with photons from the extragalactic 
background radiation fields dates baclc to the 1960s [1-4]. The reaction 77 — >■ e^e~ 
between a 7-ray of energy E and a background photon of energy e can occur when 
the center-of-mass square energy s is above threshold, s = 2Ee{l — cos 6) > 4mlc^, 
where 6 is the angle between the two photons' direction vectors. A 7-ray of energy 
i^TeV TeV therefore interacts only with background photons above a threshold en- 
ergy ethr ~ 0.3eV/i?TeV- Since the number density of background photons decreases 
roughly as a power law in energy, most of the collisions occur near threshold. Thus, 
when estimating the extinction of 1 TeV gamma rays, it is the density of the in- 
frared background which dominates; at 20 GeV, however, only UV photons near 
the Lyman limit can act as targets. 

This mechanism was recently invoked [5] to explain why many EGRET blazars 
are not seen at ~TeV energies by ground-based instruments such as Whipple, in 
spite of the fact that an extrapolation of the EGRET power-law spectra places 
them above the sensitivity limit of these ground-based detectors. The opacity r 
seen by a 7-ray in its propagation from source to Earth is roughly r ~ Naxd, 
where is the number of target soft photons above threshold, ax is the Thompson 
cross section, and d is the distance to the source. For sources with redshift z > 0.1 
(corresponding to most of the EGRET blazar sources), based on estimates of the 




FIGURE 1. Left: The mean co-moving emissivity of stellar populations as a function of red- 
shift for three different wavelengths, with (solid lines) and without (dashed lines) our metallicity 
correction added. The dashed lines are essentially a reproduction of the results of Ref. [11], and 
the observational data points arc from the Canada- French rcdshift survey group [13]. Right: The 
computed co-moving background energy density as a function of wavelength for several redshifts. 
The data points shown are high galactic latitude detections or limits at redshift z = 0: A, Ref. 
[14]; D, Ref. [15]; H, Ref. [16]; M, Ref. [17]; Te, Ref. [18]; To, Ref. [19]; V, Ref. [20]. 

diffuse IR background [6] the opacity is greater than unity for ~ TeV 7-rays, making 
their ground-based detection unlikely. 

With the advent of a new generation of ground-based instruments with antici- 
pated 7-ray energy thresholds as low as 20 GeV, and with the hkely future launch 
of GLAST [7] with sensitivity in the range ~0.01 to 100 GeV, it is important to ex- 
tend the opacity calculations down to the lowest relevant 7-ray energies. Although 
efforts along these lines have been made [8,9], very recent work on the evolution of 
star formation rates with redshift [10,11] justifies a new and more detailed calcula- 
tion of 7-ray opacity. 

The role played by the extragalactic starlight background (ESB) in the attenua- 
tion of 7-rays from extragalactic sources is defined in the exact expression for the 
7-ray opacity r, 

r{Eo,Ze) = c dz— dx- diy{l + zf^^a^^{s), 1 
Jo dz Jo 2 Jo hv 

where £'0 is the observed 7-ray energy, z^ is the source redshift, t{z) is the cosmic 
time, x = (1 — cos ^), ^ is the angle between the photons' direction vectors, v is the 

target photon frequency at redshift Uy{z) is the photon energy density per unit 
frequency at redshift z, and the Bethe-Heitler cross section for 77 ^ e+e . 



Apart from the uncertainty in cosmological parameters, the only unknown in the 
above equation is the ESB energy density Ui,{z). This can be determined if the 
mean emissivity per unit frequency, S^{z), of starhght from galaxies is known: 

u,{z)= (i^'-j-£,/(^')e"'"^'°"'^"'''^^- (2) 

Jz dz 

Here Zmax is the redshift for the turn-on of star formation, z/' = u^l+z')/ (l + z), and 
the last factor accounts for the partial absorption of the starlight by intervening 
Lya clouds during the ESB's propagation through intergalactic space [12]. 

The mean emissivity £u{z) is the total stellar energy output per unit volume 
and frequency, averaged over all galaxies and proto-galaxies, at a given redshift. 
Consider a population of stars all born at the same instant, with an initial mass 
function (IMF) (p{M) dM oc M''^ dM {a = 2.35 here). The total emission S^{T) 
from this population is the integral of the spectral energy output of each star (a 
function of its mass M, age T, and to a smaller extent its metaUicity [21,22]) 
weighted by the IMF. As the age T of the population increases, Sv{T) becomes 
redder, due to the shorter lifetimes of the bluer stars. The mean emissivity £u{z) is 
then the convolution of Su{T) with the redshift-dependent stellar formation rate, 
Ps{z): 

8,{z) = J^"^^ dz' ^Psiz')S. [T = t{z) - t{z')] C{u, z'), (3) 

where Td^g^u) is the probability that stellar photons of frequency i/ will escape 
absorption by dust and gas in their parent galaxy, and z) is a frequency- 
dependent correction to Si, which accounts for the increase in stellar metallicities 
with decrease in z. 

Figure 1 shows our results for £i,{z) (Eq.3) and u,y{z) (Eq.2). For S,y{T) we 
have used the population synthesis models of Refs. [23,24,22]; for C{iy, z) we have 
constructed an empirical correction function based on the work of Ref. [21]; the 
star formation rate Ps{z) comes from the beautiful analysis of Refs. [10,25]. (See 
Ref. [26] for more details.) 

With Eq.l, the ESB opacity to 7-rays is calculated, and shown in Fig.2, both with 
and without the metaUicity correction function £ included. Given the uncertainties 
associated with C [26], the true opacities likely lie somewhere between the two sets 
of curves. 

Figure 3 shows the effect of the ESB on 7-ray propagation from several blazars. 
Note that the spectral cutoffs occur at lower energies for blazars at higher redshifts, 
a distinctive signature which can discriminate this cutoff mechanism from intrinsic 
(intrasourcc) cutoff mechanisms. Also note that there is essentially no attenuation 
below 10 GcV, due to the sharp break in the energy density above the Lyman limit 
(Fig.l). Figure 3 also shows the beginning of the extinction of that component of 
the extragalactic 7-ray background above 20 GeV that is due to unresolved blazars 
[27] ; this is compared with recent EGRET measurements of the extragalactic 7-ray 
background [28]. 
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FIGURE 2. The opacity r of the ESB to 7-rays as a function of 7-ray energy and source redshift 
z. Left: Curves calculated with the metallicity correction included. Right: Curves calculated 
without metallicity correction. The truth likely lies between these two sets of curves. We note 
that the opacities obtained here are independent of the value assumed for Hubble's constant (see 
Ref. [26] for details). 
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FIGURE 3. Left: The attenuated power-law spectra of four prominent blazars. The solid 
(dashed) curves are calculated with (without) the metallicity correction factor. Right: Extra- 
galactic 7-ray background spectrum from unresolved blazars, calculated for the EGRET point 
source sensitivity of 10~^ cm~^s~^; solid (dashed) line includes (does not include) metallicity 
correction, and data points are from EGRET [28]. 
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